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ABSTRACT 

The  techrdque  described  presents  a  method  to  express  uncertainty 
quantitatively  in  advanced  system  cost  estimates.  In  particular,  the 
technique  suggests  the  employment  of  subjective  probability  distributions, 
which  describe  the  uncertainty  in  each  system  element,  to  determine  an 
approximate  distribution  for  total  system  cost.  A  7090  program  has  been 
written  to  perform  the  computational  operations. 
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GLOSSARY* 


Stochastic:  An  adjective  referring  to  an  operation  which  describes  uncer¬ 

tainty  statistically. 

Random 

Variable:  A  quantity  whose  value  is  uncertain.  To  fully  describe  a 

random  variable ,  the  probability  of  every  possible  value  of 
the  random  variable  must  be  given. 


Upper  (Lower) 

Tail:  The  value  of  a  random  variable  such  that  the  probability  of 

this  value  or  any  higher  (lower)  value  equals  a  stated  amount. 

Mode:  The  most  probable  value  of  a  random  variable. 


Dependency: 


If  x  and  y  are  dependent  random  variables,  knowledge  of 
the  value  of  either  variable  will  change  the  analyst's  feeling 
about  the  distribution  of  the  other.  Mathematically,  if  x  and 


y  are  dependent,  then  p  (x^/y^)  can  not  equal  p 
values  of  k  and  j. 


for  all 


Linearly  Scaled 
Beta  Distribu¬ 
tion;**  A  continuous  finite  unimodal  distribution  that  may  be  either 

skew  or  symmetrical  and  has  four  degrees  of  freedom.  The 

rectangular  distribution  is  a  special  case  of  the  scaled  Beta, 

and  the  normal  and  gamma  distributions  are  limiting  forms. 


♦Definitions  correspond  to  word  use  in  this  report. 
♦♦This  distribution  is  equivalent  to  Pearson  Type  I. 
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Mathematical  Expression: 


PBETA(x;a,M,b»=r^^1f(^)C‘(l-^)  , 


where 


.th  . 

U  Origin 

Moment: 


a,  b,  a,  and  0  are  real  numbers; 

b,  a,  and  (3  are  non -negative; 

F  (u)  is  the  gamma  function;  and 
a  <  x  <  a  +  b. 


The  average  (or  expected  value)  of  a  random  variable  raised 
th 

to  the  k  power.  The  first  origin  moment  is  the  mean ,  p. 


kth  Central 
Moment: 


Mathematical  Expression: 

(k)  V  k 

a'T  =  )  x.  p.  or 

L-,  1  1 

i 

discrete 


r+»  k 
\  X  P  (x) 

continuous 

case 


The  average  (or  expected  Value)  of  a  r;mdom  variable  dimin¬ 
ished  by  its  mean  raised  to  the  kth  power.  The  first  central 
moment  is  zero  and  the  second  is  the  variance  (the  standard 
deviation  squared). 

Mathematical  Expression: 

°r  ^"(x-^kp(x)dx  . 


discrete 

case 


continuous 

case 


Additive 

Moments;  Moments  (of  random  variables)  which  have  the  following  prop 

erty:*  the  moment  for  the  sum  of  independent  random 

til 

variables  equals  the  sum  of  the  k  moments  of  each  of  the 
added  variables. 

Mathematical  Expression  for  the  First  Four  Additive 
Moments; 


I 


Multiplicative 

Moments:  Moments  (of  random  variables)  which  have  the  following  prop¬ 

erty:*  the  k^1  moment  for  the  product  of  independent  random 
variables  equals  the  product  of  k**1  moments  Of  each  of  the 
multiplied  variables. 


Mathematical  Expression  for  Multiplicative  Moments; 

W 

Each  variable  the  analyst  defines  in  the  process  of  estimating 
the  system  total  cost.  Examples  are  quantities  Of  personnel, 
equipment  packages ,  and  computer  programs,  price  levels, 
and  planning  factors.  As  a  convenience  the  abbreviated  form, 
cost  element,  is  used  in  the  paper. 

♦ 

See  Appendix  El  for  proof. 


System  Cost 
Input  Elements: 
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Cost 

Structure: 


The  manner  in  which  input  cost  elements  are  combined  to 
determine  total  system  cost. 
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AUTHOR'S  PREFACE 


Although  there  is  wide  agreement  concerning  the  logico -mathematical 
rules  to  which  probabilities  (whatever  they  may  be)  must  adhere,  there  is 
substantial  disagreement  among  applied  mathematicians,  decision  theorists, 
and  statisticians  as  to  the  essential  meaning  of  probability.  Currently ,  this 
controversy  seems  to  be  centered' around: two  distinct  points  of  view.  One  is 
the  "objectivist"  or  "relative  frequency"  point  of  view,  which  defines  proba¬ 
bility  as  the  long-run  relative  frequency  limit  of  the  ratio  of  the  observed 
number  of  favorable  events  to  the  total  number  of  observed  instances  associ¬ 
ated  with  the  outcomes  of  a  random  physical  process.  Probability  ,  according 
to  this  point  of  view,  is  a  phenomenological  concept— a  statistic  estimated 
from  repeated  observations  of  some  directly  Observable  phenomenon.  On 
the  other  hand,  the  "subjectivist"  or  "personalistic"  point  of  view  defines 
probability  as  a  numerical  coefficient  purporting  to  measure  a  particular 
human  being's  subjective  belief  about  the  outcomes  of  some  physical  phenom¬ 
enon.  Probability ,  from  the  subjectivist  standpoint,  is  not  a  phenomenologi¬ 
cal  concept  at  all.  It  is  not  a  characteristic  of  the  physical  phenomenon  to 
which  it  purports  to  refer,  but,  rather  ,  a  characteristic  of  human  beings— a 
component  part  of  a  particular  individual's  attitude  toward  a  physical  phe¬ 
nomenon.  As  such,  subjective  probability  is  a  fictional  concept,  much  like 
the  concept  of  "force"  in  physics,  which  can  only  be  inferred  from  observa¬ 
tions  of  displacement  and  motion  on  the  part  of  physical  bodies,  and  like  the 
concept  of  "intelligence"  in  psychology,  which  can  only  be  inferred  from 
observed  verbal  and/or  choice  behavior  on  the  part  of  human  beings. 

Viewing  probability  from  the  objectivist  point  of  view  has  proved  highly 
successful  in  many  historical  applications,  including  the  management  of 
gambling  houses,  industrial  quality  control,  and  numerous  scientific  endeav¬ 
ors,  particularly  in  the  field  of  genetics.  However,  the  success  of  these 
applications  rests  upon  the  existence  of  a  stable,  physical  process  of  which 
repeated  observations  can  be  made.  Even  more,  the  whole  notion  of  objec¬ 
tive  probability  requires  the  existence  of  such  processes,  since  the  objectiv¬ 
ist  definition  does  not  apply  to  a  nonrepetitive  phenomenon.  Strictly  inter¬ 
preted,  therefore,  the  objectivist  point  of  view  can  be  of  no  assistance  in 
explaining  or  predicting  the  outcome  (s)  of  nonrepetitive  phenomena. 

However,  many  real-world  decision  makers  face  the  necessity  of 
taking  definite  actions  in  the  face  of  substantial  uncertainty  regarding  the 


outcomes  of  nonrepetitive  phenomena.  Since  the  objectivist  concept  of  proba¬ 
bility  cannot  be  of  service  in  making  such  decisions,  the  subjectivist  view 
has  been  developed  to  fill  the  void.  The  conceptual  foundations  of  subjective 
or  personalistic  probability  have  been  worked  out  by  Leonard  J.  Savage.  M* 
Two  pioneers  in  applying  this  concept  to  many  real-world  decisions  problems 
have  been  Robert  SchlaifferC2!  and  Howard  Raifia,  [3]  Specific  implementa¬ 
tion  of  this  concept  in  oil  well  drilling  decisions  is  reported  by  C.  Jackson 
Grayson. 

The  position  taken  in- this  report  is  that  the  process  by.  which  military 
system  costs  are  generated  does  not  usually  constitute  a  stable,  repetitive 
phenomenon  to  which  objectivist  probability  is  an  applicable  concept. 
Therefore,  probability  is  considered  from  the  subjectivist  point  of  view, 
and  the  author  presents  a  technique  by  which  the  necessary  implications  of 
a  decision-maker’s  or  analyst's  subjective  judgments  concerning  costs  may 
be  derived  and  displayed. 


J.  R.  Miller,  ffl 


♦Numbers  in  brackets  refer  to  references  at  end  of  report. 


SECTION  I 


INTRODUCTION 


Ite, 


PURPOSE 

The  analyst  faced  with  the  task  of  estimating  the  cost  of  advanced  mili¬ 
tary  systems  is  continually  beset  by  the  problem  of  coping  with  a  host  of 
attendant  uncertainties.  Single -value  estimates,  although  duly  qualified,  belie 
the  full  range  of  possible  cost  consequences  that  may  ensue  once  the  decision 
to  acquire  a  given  system  has  been  made,  and  in  a  great  many  instances  such 
estimates  gravely  misstate  the  ultimate  system  cost.  Therefore,  early  in 
the  system  life  cycle ,  defense  planners  have  been  promulgating  their  need  for 
improved  information  about  a  proposed  system's  possible  economic  implica¬ 
tions,  despite  manifest  deficiencies  in  any  early  system  description.  The 
present  paper  suggests;  the  kinds  of  information  about  the  uncertainty  which 
might  be  shown  to  defense  planners  and  offers  a  technique  to  generate  such 
information. 


SOURCES  GF  UNCERTAINTY 

Uncertainties  in  advanced  systems  cost  estimates  arise  from  many 
Sources.  Examination  of  such  sources  has  revealed  that  they  may  be  divided 
into  two  basic  categories.  The  first  category,  examples  of  which  are  shown 
in  Table  I,  consists  of  those  types  of  uncertainty  which  are  amenable  to  plan¬ 
ning  decisions,  that  is,  those  types  over  which  decision  makers  exercise 
control.  The  second  category,  examples  of  which  are  shown  in  Table  n,  is 
the  class  of  uncertainties  which  is  beyond  direct  control  or  explicit  planning 
during  the  conceptual  phase  of  a  system. 
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EXPRESSING  COST  UNCERTAINTY 


Because  of  fundamental  differences  between  the  two  categories  of  uncer¬ 
tainty,  different  techniques  are  needed  in  dealing  with  them.  Category  I 
uncertainties  might  best  be  represented  by  exhaustively  identifying  ail  signifi¬ 
cant  competing  design  alternatives  and  by  providing  the  decision  maker  with 
the  economic  implications  of  each  planning  choice.  Nevertheless ,  because  of 
the  large  lead  time  between  the  generation  of  advanced  system  cost  estimates 
and  the  date  by  which  the  system  achieves  operational  status,  a  large  number 
of  residual  uncertainties,  i.e. ,  Category  n  uncertainties,  will,  in  general, 
characterize  each  of  the  identified  system  alternatives.  It  would  thus  be 
desirable  to  convey  to  the  decision  maker  a  meaningful  numerical  range 
describing  the  magnitude  of  the  effect  on  the  total  system  cost  of  the  simul¬ 
taneous  presence  of  all  the  residual  uncertainties  in  the  system  elements.  In 
addition,  numerical  ranges  describing  lower  levels  of  Cost  aggregations  than 
the  total  system  are  valuable  to  expose  major  sources  of  uncertainty .  Finally , 
to  support  statistical  analysis,  it  would  be  useful  to  show  the  relative  likeli¬ 
hood  that  different  cost  magnitudes  will  be  the  actual  ultimate  total  system 
cost.  This  report  is  addressed  to  generating  this  information. 

Table  I 

Examples  of  Sources  of  Category  I  Uncertainty 
System  Performance  Requirements  * 

Examples:  coverage,  numbers  of  users  served,  survivability,  reliability, 
flexibility,  growth  potential,  reaction  time ,  range ,  warning  time. 

_  C  . 

♦Although  these  requirements  depend  on  factors  (such  as  the  enemy  threat  and 
existing  technology)  which  a  decision  maker  does  not  control,  he  does  generate 
an  explicit  and  conscious  statement  of  the  performance  requirements  for  the 
system  on  the  basis  of  a  given  set  of  constraints  and  goals;  thus  system  per¬ 
formance  requirements  are  included  in  Category  I  uncertainties. 
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Implementation  Plan 

Examples:  schedules,  amount  of  in-house  vs.  external  management, 
number  and  extent  of  competing  contractor  studies  during  the  program  defi¬ 
nition  phase. 

Configuration  (Resource  inventory) 

Examples:  number  and  kinds  of  hardware  packages,  the  nature  and  extent 
of  software,  and  the  quantity  and  quality  of  personnel.  ** 

Operational  Concept 

Examples:  modes  of  operation  during  varying  periods  (e.  g. ,  tranquil, 
crisis  ,  tfans  and  post  attack) ,  amount  of  peace  time  exercises. 

Maintenance  Concept 

Examples:  contractor  vs.  military  maintenance ,  number  of  shifts. 
Logistics- Concept 

Examples:  number  ,  size,  and  locations  of  warehouses  to  maintain  supply 
line ,  quantities  and-types  of  stocks  and  spares  to  be  stored  at  each; 

Training  Concept 

Examples:  in-house  vs.  contractor  training,  formal  vs.  on-the-job 
training. 

Funding  Schedule 

Examples:  scheduling  of  obligational  authority. 


♦•Ordinarily  in  the  conceptual  phase ,  these  can  be  only  grossly  defined. 
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Table  n 


Examples  of  Sources  of  Category  II  Uncertainty 
Price  levels  for  equipment  and  contractor  services 

Extent  of  R&D  effort  needed  for  desired  equipment  development  and  resulting 
production  Cost  of  developed  items 

Exact  quantities  and  types  of  personnel  required  to  effectuate  a  given  opera¬ 
tions  and  maintenance  concept 

Exact  quantities  and  types  of  support  equipment  (AGE) ,  initial  and  follow-on 
spares,  documentation,  etc. 

Manpower  and  computer  time  requirements  for  system  design  and  management 

Exact  hardware  and  software  design  specifications  for  operational  system 

A  FEW  EARLY  APPROACHES  TO  DESCRIBE  CATEGORY  II  UNCERTAINTIES 

A  common  practice  in  cost  estimating  is  to  study  the  nature  and  magni¬ 
tude  of  the  uncertainty  in  each  of  the  system  cost  input  elements*  (See 
Glossary)  and  footnote  the  estimate  with  statements  about  possible  variations 
imtheir  values.  This  procedure  serves  to  put  the  user  of  a  cost  estimate  on 
notice  that  the  single  value  shown  for  total  system  cost  is  Subject  to  error, 
and  helps  to  focus  attention  bn  the  major  sources  of  uncertainty.  Neverthe¬ 
less,  such  information  does  not  reveal  the  extent  to  which  the  estimated  total 
system  cost  is  likely  to  deviate  from  the  actual  cost. 

In  order  to  provide  some  estimate  of  uncertainty  in  total  system  cost, 
analysts  have  attempted  to  make  an  astute  guess  about  the  variability  of  Such 

♦For  simplicity ,  the  abbreviated  form,  cost  element,  will  be  used  subse¬ 
quently  throughout  this  report. 


cost  and  record  the  result  for  the  decision  maker’s  reference.  The  analyst 
performs  this  operation  by  mentally  assessing  the  combined  impact  of  all 
sources  of  uncertainty  in  all  of  the  cost  elements.  For  a  typical  size  system, 
such  practice  strains  the  faculties  of  even  the  most  gifted  analyst,  and  hence 
the  results  are  of  questionable  validity. 

Another  expedient  used  in  quest  of  more  reliable  information  about  the 
variability  in  total  costs  requires  the  analyst  to  specify  a  likely  range  for  each 
of  the  cost  elements;  then  two  additional  total  system  cost  estimates  can  be 
computed,  one  being  the  sum  of  the  lowest  values  for  all  the  cost  elements , 
and  the  other  being  the  sum  of  the  highest  values.  In  this  manner  a  range  iS 
established  for  the  total  System  cost.  However,  such  a  range  is  subject  to 
strong  criticism.  It  is  Self-evident  that  the  possibility  that  all  elements  will 
actually  attain  their  lowest  values  (or  their  highest)  is  very  remote.  Hence, 
Such  a  range  is  a  serious  overstatement  of  the  magnitude  of  likely  variability 
in  total  system  costs. 

A  more  precise  technique  is  therefore  needed  to  translate  the  analyst's 
feelings  about  the  uncertainty  in  all  the  cost  elements  into  a  statement  about 
the  uncertainty  in  the  total  system  cost.  This  paper  has  been  written  to 
suggest  one  such  technique. 


SECTION  II 


DESCRIPTION  OF  THE  PROPOSED  TECHNIQUE 


GENERAL 

Uncertainty  in  each  of  the  cost  elements  is  described  by  treating  the 

C2] 

elements  as  random1  variables.  A  subjective  distribution '  J  for  each  element 
is  provided  by  the  analyst  to  represent  his  feelings  about  the  relative  likeli¬ 
hood  of  all  feasible  values.  The  burden  placed  on  the  analyst,  however,  is 
minimized  by  requesting  that  he  provide  only  four  numbers  which  are  esti¬ 
mates  of  four  parameters  characterizing  the  shape  of  the  distribution  curve. 


From  the  four  parameters  the  entire  distribution  is  inferred.  This  is 
done  by  assuming  a  linearly  scaled  beta  function  (see  Glossary)  which  satis¬ 
fies  the  analyst's  inputs.  The  element  distributions  are  then  combined  in 

[5]  r 

accordance  with  a  cost  structure  L  to  produce  the  probability  distribution 


for  the  total  system  cost  and  equivalent  distributions  for  intermediate  levels 


of  cost  aggregation.  Cumulative  probability  intervals  and  relevant  summary 


statistics  (e.g. ,  mode,  mean,  or  median)  can  be  generated  from  such  dis¬ 


tributions. 


LIMITATIONS 

The  procedure  described  is  suited  for  direct  application  in  the  following 
situation:  (1)  the  probability  distributions  specified  in  the  input  data  are  uni- 
modal;  (2)  only  additive  and  multiplicative  operations  are  involved  in  relating 
element  costs  to  total  costs;  (3)  elements  which  are  multiplied  together  are 
independent;  and  (4)  elements  which  are  added  together  are  either  independent 
or  linearly  dependent  (see  Appendix  1). 


Procedure  are  presented  for  permitting  application  of  the  technique 
in  some  situations  which  do  not  satisfy  the  above  requirements. 

INPUTS  PROVIDED  BY  THE  COST  ANALYST 

frl 

To  cope  with  uncertainty  in  schedule  estimates,  PERT1  J  has  adopted 
a  procedure  which  utilizes  a  scaled  beta  function  and  requires  that  the  analyst 
furnish  three  numbers  to  specify  its  form.  These  three  numbers  are  esti¬ 
mates  of  the  lowest,  X  ,  and  highest,  X  ,  values  (lower  and  upper 

L  H 

1  percent  tails,  respectively)  and  the  mode,  Xp  ,  (most  probable  value). 

The  PERT  procedure  assumes  that  the  distribution’s  standard  deviation 

can  always  be  satisfactorily  approximated  as  one-sixth  of  the  total  range, 

i.e. ,  1/6  X  -X  .  However,  since  1  percent  tails  are  difficult  to  esti- 
n  L 

mate  accurately,  standard  deviations  calculated  in  this  manner  may  be  unre¬ 
liable;  in  addition,  this  method  of  computation  makes  no  allowances  for 
widely  dispersed  or  sharply  peaked  distributions.*  Hence,  the  resulting 
curve  often  grossly  distorts  the  analyst's  actual  feelings.  Furthermore,  this 
restriction  on  the  standard  deviation  seriously  reduces  the  variety  of  shapes 
which  are  inherent  in  the  beta  family. 

To  overcome  such  shortcomings ,  it  has  been  decided  in  the  proposed 

technique  to  specify  a  fourth  number,  the  80  percent  central  range.  The 

80  percent  central  range ,  C  ,  represents  the  range-  of  the  estimated 

R 

variable  which  contains  80  percent  of  the  probability.  It  is  calculated  by 

subtracting  the  lower  end  of  the  range  (the  lower  10  percent  tail)  from  the 

upper  end  of  the  range  (the  upper  10  percent  tail).  By  permitting  the  analyst 

the  freedom  of  specifying  C  ,  more  representative  distributions  often  can 

R 


♦Statisticians  refer  to  this  phenomenon  as  the  distribution’s  kurtosis. 


Figure  1 .  Symmetrical  Beta  Distributions 
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Figure  2.  Skew  Beta  Distributions 
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result.  For  example ,  the  rectangular ,  normal ,  and  gamma  distributions  may 

* 

all  be  accurately  portrayed.  Figures  1  and  2  are  provided  to  demonstrate 
some  of  the  varied  shapes  that  t!;e  beta  distribution  can  assume. 

Hence,  for  each  cost  element,  the  analyst  provides  four  values,  XT  , 

**  L 

Xjj  ,  Xp  ,  and  ,  to  describe  the  uncertainty.  This  information  is 

based  on  the  analyst's  assessment  o£  the  probable  sources  and  magnitude  of 
the  uncertainties.  For  such  assessment,  the  analyst  not  only  uses  relevant, 
processed,  historical  data,  but  obtains  and  evaluate*  the  opinions  of  experts 
and  takes  into  account  any  discrepancies. 

At  this  point  it  might  be  well  to  discuss  the  analyst's  willingness  to 
provide  more  data  than  in  the  past.  It  is  quite  probable  that  the  analyst's 
reaction  to  such  a  situation  will  be  related  to  his  integrity.  If  he  sincerely 
believes  that  extensive  uncertainties  enshroud  the  value  of  a  particular  cost 
element,  it  will  be  very  difficult  for  him  to  specify  a  single  value.  Hence, 
it  is  the  author's  contention  that  when  uncertainty  exists,  the  conscientious 
cost  analyst  will  appreciate  the  opportunity  of  being  able  to  stipulate  ranges. 

DETERMINING  UNCERTAINTY  IN  TOTAL  COSTS 

Random  variables  representing  cost  elements  may  be  combined  in 
accordance  with  a  cost  structure  to  produce  a  random  variable  to  represent 
the  total  system  cost.  By  knowing  (1)  the  distributions  of  the  cost  elements, 

*  The  linearly  scaled  beta  function  is  perhaps  more  generally  known  as  a 
Pearson  I  curve.  The  gamma  and  normal  are  limiting  forms  and  hence  can 
be  approximated  with  any  desired  accuracy.  The  rectangular  distribution  is 
a  special  case  of  the  beta. 

**When  the  uncertainty  is  insignificant,  only  one  value  need  be  specified. 
This  value  should  be  set  equal  to  the  ro.vva  and  all  higher  moments  can  be 
taken  equal  to  zero. 
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(2)  the  nature  of  any  statistical  dependence  which  may  exist  among  them,  and 

(3)  the  functional  relationships  which  reflect  the  cost  structure ,  the  distribu¬ 
tion  for  the  total  system  cost  may  be  determined. 

To  satisfactorily  treat  uncertainty  in  advanced  system  costs,  it  is  con¬ 
sidered  adequate  to  represent  the  cost  structure  by  products  and  sums  of 
random  variables  where  multiplied  variables  are  independent  and  added  vari¬ 
ables  are  either  independent  or  linearly  dependent.  The  concept  of  linear 
dependence,  as  used  in  this  paper,  is  defined  in  Appendix  I. 

The  additive -multiplicative  restriction  on  the  cost  structure  is  not 
believed  to  be  a  very  serious  handicap,  since,  in  the  author’s  experience , 
these  operations  are  most  frequently  encountered  in  costing  advanced  elec¬ 
tronic  command  and  control  systems.  However,  it  is  possible  to  process 
more  general  functional  relationships,  as  described  in  Appendix  II. 

The  limitations  placed  on  intervariable  dependence  are  also  believed 
not  to  be  highly  restrictive.  Appendix  I  shows  that  in  many  situations  depend¬ 
ence  can  be  avoided  completely  by  the  proper  choice  of  cost  element  vari¬ 
ables.  However,  in  those  instances  where  this  expedient  is  not  possible, 
linear  dependence,  in  the  author’s  opinion,  will  generally  be  adequate  to 
reflect  the  statistical  relationship  among  variables. 

As  demonstrated  in  Appendix  I,  the  sum  of  linearly  dependent  variables 
may  be  expressed  as  sums  and  products  of  independent  auxiliary  random 
variables.  The  appendix  further  describes  the  information  which  the  analyst 
may  provide  to  determine  these  auxiliary  variables.  Thus,  in  a  great  many 
circumstances,  the  author  believes  that  the  cost  structure  may  be  satisfac¬ 
torily  represented  by  the  sums  and  products  of  independent  random  variables. 
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It  is  well  known  that  when  independent  random  variables  are  added  the 
distribution  of  their  sum  can  be  determined  by  convolving  the  component 
distributions.  Similarly,  when  independent  random  variables  are  multiplied, 
the  distribution  of  their  product  can  be  determined  by  an  analogous  process 
(which  is  portrayed  in  Appendix  n) .  Nevertheless,  since  the  input  data  about 
the  uncertainty  are  only  approximate,  it  was  felt  that  such  elaborate  opera¬ 
tions  were  not  warranted.  Instead,  the  decision  was  made  to  use  moments 
to  summarize  the  distributions.  This  procedure  is  discussed  in  the  following 
section. 

UTILIZATION  OF  MOMENTS  TO  GENERATE  THE  DISTRIBUTION  OF  TOTAL 
COSTS 

To  treat  uncertainty  in  PERT,  the  mean  and  standard  deviation  for  each 
element  distribution  is  computed  and  then  used  to  compute  the  mean  and 
standard  deviation  for  the  total.  The  central  limit  theorem,  which  states 
that  Under  suitable  conditions  the  distribution  representing  the  sum  of  random 
variables  tends  to  the  normal,  is  then  customarily  invoked.  Since  the  normal 
distribution  is  completely  specified  by  the  mean  and  standard  deviation,  the 
distribution  of  the  total  is  thus  determined. 

The  degree  to  which  the  normal  faithfully  portrays  the  distribution  of 
the  sum  depends  on  many  factors ,  including  the  number  of  variables  added, 
the  shapes  and  relative  sizes  of  their  distributions ,  and  the  degree  of 
dependence  among  them.  It  is  not  clear,  therefore,  that  total  system  costs 
will  always  be  adequately  represented  by  the  normal,  and  it  is  even  less 
apparent  that  the  distribution  of  intermediate  levels  of  cost  aggregation  will 
be  adequately  represented.  The  mean  and  standard  deviations  themselves 
provide  no  clue  on  the  reliability  of  the  normal  approximation. 
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Figure  3  provides  a  graphic  portrayal  of  the  kinds  of  probability  state¬ 
ments  that  can  be  made  for  a  normal  distribution.  It  also  shows  the  upper 
bound  for  confidence  statements  that  can  be  made  for  any  distribution.  In 
addition,  it  depicts  an  estimated  upper  bound  for  unimodal  distributions. 

Because  confidence  statements  for  the  normal  are  so  much  tighter  than 
for  either  of  the  indicated  upper  bounds ,  use  of  the  latter  bounds  will  fre¬ 
quently  grossly  overstate  the  uncertainty.  Furthermore ,  the  mean  and 
standard  deviation  provide  no  insight  about  skewness.  Thus ,  it  is  desirable 
to  generate  more  information  about  the  uncertainty  in  cost  aggregations 
before  attempting  to  construct  a  distribution. 

To  cope  with  the  above  situation,  it  was  decided  to  use  the  first  four, 
instead  of  just  the  first  two,  moments  (such  as  the  mean  and  standard  devi¬ 
ation).  In  this  regard,  additive  and  multiplicative  moment  sets  were 
employed.  These  are  defined  in  the  glossary  and  their  properties  are 
demonstrated  in  Appendix  III.  The  n^  idditive  moment  for  a  sum  of  inde¬ 
pendent  random  variables  is  simply  computed  by  adding  the  n*1  additive 
moments  for  each  of  the  variables.  Equivalently,  the  n^  multiplicative 
moment  for  a  product  of  independent  random  variables  is  the  product  of  the 
n**1  multiplicative  moments  of  each  variable. 

The  distribution  for  total  system  cost  (and  for  intermediate  levels  of 
cost  aggregation)  is  constructed  by  determining  the  linearly  scaled  beta* 

(a  four -parameter  function)  which  has  the  same  first  four  moments.  The 
various  mathematical  operations  that  are  involved  in  this  process  are 


*If  the  beta  distribution  cannot  satisfy  these  moments ,  the  appropriate 
member  function  from  the  Pearson  family  may  be  used  (see  Table  43, 
Ref.  [7]). 
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Figure  3.  Probability  that  a  Random  Variable  is  More  than  Ko  Away 
from  the  Mean 
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described  in  Appendix  IV.  The  integrity  of  this  procedure  is  demonstrated  in 
Appendix  V  by  comparing  the  resultant  distributions  with  the  theoretical 
curves. 
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SECTION  m 
COMPUTER  PROGRAM 

*■ 

r 

j  INTRODUCTION 

i 

(  A  program  has  been  written  for  the  7090  to  facilitate  the  translation  of 

-  the  analyst's  feelings  about  the  uncertainty  in  each  cost  element  into  the 

;  parameters  of  an  approximating  linearly -scaled  beta  distribution  for  the  cost 

total  and  sub-totals.  These  parameters  are  used  to  generate  upper  and 
lower  bounds  to  a  cumulative  probability  interval  of  any  desired  size.  The 
program  is  composed  of  three  discrete  parts  which  are  described  below. 

;  The  reader  is  referred  to  Appendix  VI  for  a  more  detailed  description 

|  of  the  program,  including  an  explanation  of  the  mechanics  involved  to  operate 

;  it  on  the  MITRE  7090  Facility. 

{ 

PROGRAM  I 

Program  I  accepts  as  inputs  those  four  parameters  of  the  linearly 
j  scaled  beta  distribution  which  are  provided  by  the  analyst.  The  reader  will 

recall  that  these  are  the  most  probable  value,  the  lowest  and  highest  values, 
and  the  80  percent  central  range.  With  this  information,  the  four  additive 
moments  are  computed  and  presented  as  an  output.  The  input  is  put  on 
punch  cards  and  the  output  is  printed.  The  reader  is  referred  to  Figure  4 
for  a  graphic  description  of  this  program. 

PROGRAM  H 

Program  Q  accepts  the  additive  moments  computed  by  Program  I  for 
all  of  the  cost  elements  and  computes  an  equivalent  set  of  moments  for 
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LINEARLY  SCALED 

x,  y,  v,  a  w 

SUBJECTIVE  ESTIMATES 
PROVIDED  BY  ANALYST 

BETA  DISTRIBUTION 

(FIRST  FOUR 
MOMENTS) 

ADDITIVE 

FOR  EACH  COST  ELEMENT 


X,  MOST  LIKELY  VALUE  (MODE) 
xL  a  xH  LOWER  a  UPPER  i%  tails 

C„  80%  CENTRAL  RANSE  (UPPER  10%  TAIL  MINUS  LOWER  10%  TAIL) 

X  *  1st  ORIGIN  MOMENT  *  MEAN  (EXPECTED  VALUE) 

Y  *  2nd  CENTRAL  MOMENT  =  VARIANCE  (STD  DEVIATION*) 

V  =  3rd  CENTRAL  MOMENT 

W  s  4th  CENTRAL  MOMENT  minut  Y* 

Figure  4.  Graphic  Representation  of  Program  1 

l 

cumulative  costs  at  all  desired  levels  of  aggregation.  An  effort  was  made  to 
keep  the  program  sufficiently  flexible  so  that  a  large  variety  of  arithmetic 
arrangements  in  the  cost  structure  could  be  handled.  The  input  is  inserted 
on  punch  cards.  The  output,  specifying  the  additive  moments  of  aggregated 
costs,  is  furnished  on  punch  cards  (which  can  be  used  as  an  input  to  Program 
IU)  and  is  also  printed.  Figure  5  depicts  the  program. 

PROGRAM  m 

Program  in  converts  the  additive  moments  computed  by  Program  II 
into  the  parameters  (see  Glossary)  for  the  linearly  scaled  beta  distribution. 
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x,  y.  v,  a  w 

- - - ► 

FOR  EACH  SUBSYSTEM  AND 

FOR  TOTAL  SYSTEM  COST 


[RESTRICTION  ON  COST  STRUCTURE: 

SUMS  AND  PRODUCTS  OF  INDEPENDENT  RANDOM  VARIABLES] 
Figure  5.  Graphic  Representation  of  Program  11 

It  further  determines  bounds  for  any  y  cumulative  probability  interval  by 
finding  limits  for  the  1-  -  y  tails.*  Inputs  are  on  punch  cards  and  the 
output  is  printed.  Figure  6  depicts  this  program. 


♦Exact  upper  and  lower  bounds  to  these  tails  can  be  computed  directly  from 
the  first  four  moments  to  reveal  the  maximum  error  in  the  beta  approxima¬ 
tion.  [8] 
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SECTION  IV 


USEFULNESS  OF  THE  RESULTS 


Information  provided  by  the  technique  proposed  herein  furnishes  a  better 
single  value  (as  compared  to  conventional  estimates)  to  represent  the  most 
likely  cost  of  a  military  system  and,  in  addition,  furnishes  a  quantitative 
statement  about  the  nature  and  extent  of  the  uncertainty. 

Conventional  estimates  typically  introduce  the  most  likely  value  for 
each  cost  element  into  the  cost  structure  to  determine  the  system's  total 
economic  implications.  Nevertheless,  since  the  distributions  for  the  cost 
elements  are  generally  skewed  upwards ,  the  single  value  thus  computed  for 
the  system  often  has  very  little  likelihood  of  occurring  and  may  constitute  a 
serious  understatement  of  the  true  cost.  When  a  single  value  is  required 
from  the  cost  analyst,  the  mean  or  the  modal  value  (which  are  frequently 

[2] 

relatively  close  for  the  total  system  cost)  is  generally  a  more  meaningful 
figure  to  specify.  In  systems  which  are  characterized  by  considerable 
research  and  development,  these  latter  values  have  been  shown  to  exceed  the 
conventional  estimate  by  more  than  30  percent. 

Bracketed  probability  intervals  provide  a  quantitative  measure  of  the 
range  of  the  uncertainty  in  aggregate  costs.  Such  intervals  provide  insight  on 
whether  or  not  the  accuracy  of  a  cost  estimate  satisfies  the  requirements  of 
its  intended  application.  The  decision  maker  may  establish  in  advance  the 
accuracy  criterion  which  must  be  met.  The  bracketed  interval  would 
disclose  whether  the  criterion  is  satisfied. 


Bracketed  probability  intervals  for  subsystem  costs  and  other  inter¬ 
mediate  levels  of  aggregation  serve  to  identify  areas  of  major  cost  uncer¬ 
tainty.  These  are  areas  which  might  benefit  from  subsequent  efforts  to 
refine  costs.  However,  consideration  should  be  given  to  the  sensitivity  of 
the  uncertainties  to  more  intensive  study,  since  substantial  additional 
probing  will  not  materially  reduce  certain  types  of  uncertainties. 

The  probability  function  can  be  used  to  make  computations  which  may 

provide  extremely  useful  information  to  the  decision  maker.  For  example , 

when  two  or  more  equal  effectiveness,  alternative  configurations  are  being 

considered  for  acquisition,  one  can  either  immediately  select  the  system 

with  the  lowest  expected  cost  or  defer  a  decision,  pending  the  collection  of 

more  information.  By  knowing  the  entire  probability  function,  it  is  possible 

to  compute  (1)  the  probability  that  the  alternative  having  the  lowest  ultimate 

cost  will  not  be  selected  by  using  the  lowest  mean  criterion,  and  (2)  the 

expected  size  of  the  cost  premium  which  results  from  not  having  perfect 

[2] 

information.  This  expected  loss  (premium)  can  be  used  as  an  index  of 
the  maximum  effort  that  should  be  devoted  to  refining  cost  or  to  pursuing 

[9] 

exploratory  R&D.  Conversely,  this  expected  loss  represents  the  maxi¬ 
mum  expected  gain  from  deferring  immediate  action. 
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SECTION  V 

!  VALIDITY  OF  RESULTS 

j  The  reader  is  no  doubt  well  aware  of  the  fact  that  a  conventional  cost 

<  estimate  is  no  better  than  the  quality  of  its  inputs.  The  analyst's  ability  to 

|  identify  all  cost  contributing  elements  and  to  properly  assess  their  economic 

t 

|  implications  bears  directly  on  the  validity  of  the  end  product.  Similarly,  in 

!  appraising  the  sources  and  magnitudes  of  the  uncertainties  which  influence 

» 

I  each  element's  cost,  there  is  no  satisfactory  substitute  for  reliable  data. 

>  The  cost  analyst  must  carefully  and  astutely  search  out  all  significant  factors 

!  which  give  rise  to  system  element  cost  variability  and  reflect  the  magnitude 

|  of  their  effect.  The  author  wishes  to  point  out,  however,  that  although  this 

does  place  an  added  burden  on  the  analyst,  he  will  probably  derive  greater 
satisfaction  from  pursuing  this  task  than  trying  to  generate  a  single  number 
in  the  face  of  imperfect  information. 

It  is  true  that  the  analyst's  inputs  to  describe  the  uncertainty  are  sub¬ 
jective.  Nevertheless ,  because  of  the  general  unavailability  of  completely 
relevant  historical  data  from  which  the  cost  of  advanced  systems  can  be 
accurately  gauged,  single  value  estimates,  themselves,  contain  varying 
degrees  of  subjectivity.  What  is  perhaps  more  important,  however,  is  the 
fact  that  although  cost  analysts  in  the  past  have  been  fully  aware  of  the 
presence  and  character  of  element  uncertainties ,  the  means  did  not  exist  to 
properly  assess  their  cumulative  impact  on  the  total  system  cost.  Hence , 
more  often  than  not,  only  an  extremely  unreliable  conjecture  about  total  cost 
uncertainty  could  be  passed  on  to  decision  makers.  By  applying  the  proposed 
1  technique,  the  cost  analyst  is  permitted  to  focus  his  attention  on  appraising 
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the  uncertainties  surrounding  the  individual  cost  elements ,  an  area  where  he 
is  believed  more  qualified  to  judge  uncertainty;  and  standard  statistical 
methods  are  then  employed  to  measure  their  aggregate  implications  on  the 
system  total. 

If  the  analyst  inadvertently  neglects  to  consider  significant,  sources  of 
uncertainty  or  has  even  forgotten  to  include  some  cost  elements,  the  results, 
naturally,  would  understate  the  magnitude  of  the  uncertainty  and  bias  the  dis¬ 
tribution.  Oversight  of  relevant  cost  areas  is ,  of  course ,  an  error  which 
has  just  as  serious  an  effect  on  single-value  estimates,  and  hence,  much 
work  has  been  and  is  being  done  to  minimize  this  possibility.  Similarly, 
work  is  needed  in  developing  a  capability  to  realistically  appraise  the  sources 
and  magnitude  of  uncertainties. 

It  is  the  author  's  contention  that  any  information  the  analyst  can  provide 
regarding  the  magnitude  of  the  uncertainty  (even  if  occasionally  in  the  form  of 
an  educated  guessj  is  better  than' no  knowledge  at  all.  The  proposed  technique, 
in  its  present  form,  is  regarded  as  merely  a  first  step  in  the  direction  of 
providing  decision  makers  with  more  meaningful  information  about  the  vari¬ 
ability  in  total  system  costs. 
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SECTION  VI 

AREAS  OF  FUTURE  WORK 

Many  refinements  of  the  technique  described  herein  will  probably 
germinate  from  the  experience  gained  in  its  successive  application.  Two 
specific  areas  where  more  intensive  study  seems  to  offer  large  potential 
rewards  are:  (1)  development  of  techniques  to  aid  in  obtaining  quantitative 
information  about  the  cost  uncertainty  from  experts;  and  (2)  development  of 
analytical  tools  that  will  permit  fuller  utilization  of  the  generated  statements 
about  the  uncertainty  in  making  planning  decisions. 


'**}*•.. 
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APPENDIX  I 

DEPENDENCE  AMONG  VARIABLES 
AVOIDANCE  OF  DEPENDENCY 

In  a  typical  cost  estimate,  there  are  certain  cost  items  which  one 
normally  would  expect  to  be  dependent.  An  example  is  mission  equipment 
cost,  initial  spares  cost,  and  maintenance  cost.  Nevertheless,  if  the  total 
of  these  three  items  is  determined  as  below 

NxDc[1+.°1Pm(1+%)]  , 

where 

N  =  number  of  units; 

Uc  =  equipment  cost; 

PM  =  percentaSe  maintenance  costs  are  of  equipment  costs;  and 
Yg  =  ratio  of  spare  costs  to  maintenance  costs, 

the  component  variables  in  the  above  expression  may  frequently  be  regarded 
as  independent.  When  the  above  technique  can  be  employed,  it  is  possible  to 
avoid  other  means  of  specifying  the  nature  of  the  dependency. 

LINEAR  DEPENDENCE* 

Definition  and  Properties 

When  two  cost  elements  are  dependent,  the  relationship  between  the 
elements  can  often  be  satisfactorily  expressed  by  the  following  equality: 


*  The  reader  should  not  confuse  this  concept  with  die  term  as  applied  to  sets 
of  matrices  or  polynomials.  Dependence  here  is  concerned  with  the  statis¬ 
tical  relationships  among  random  variables;  see  Glossary. 
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s-(V*)( xr>Ji+Kb)+Kc 


where 


and  x.. 


Hi 

Ka*  Kb-  Kc 


=  random  variables  representing  the  cost 
elements; 

=  a  random  variable  of  zero  mean  and 
independent  of  ; 

-  the  mean  of  x1  ;  and 
=  constants. 


In  the  above  situation,  x  and  x  are  said  to  be  (for  the  purposes  of  this 

i-  u 

report)  linearly  dependent.  The  conditional  distribution  for  x  given  x  , 

A  X 

has  the  following  properties:  (1)  when  x1  equals  its  mean,  the  distribution 
for  x  can  be  specified  without  constraint;  and  (2)  the  distributions  for  other 

it 

values  of  are  represented  by  a  linear  shift  and  expansion  (or  contraction) 
of  the  coordinates.*  The  reader  will  note  that  the  sum  of  linearly  dependent 
random  variables  x  and  x  becomes: 

VV(V1  +  0(VVmi>  VVmi  ' 

which  is  composed  of  independent  random  variables. 

Determination  of  Independent  Variables 

The  variables  K&  ,  ,  Kc  and  e  are  determined  on  the  basis  of 

information  furnished  by  the  analyst.  Although  there  is  more  than  one  maimer 


♦Thus 


where 


p[VCv-0]-'<S 


b+l  /  1 


~  Ka  (xi  '  fa  '  i^  (X1  "  0 
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of  eliciting  the  necessary  information,  only  one  will  be  considered  in  this 
report. 


The  analyst  is  called  on  to  describe  the  conditional  distribution  for 


Xg  .  given  that  x^  is  at  its  mean  value,  ^  .  He  can  do  this  by  specifying 
values  for  the  parameters  of  the  linearly  scaled  four-parameter  beta.  From 
these  parameters,  the  four  additive  moments  A  ,  A  ,  A  and  A  can 

i  a  u  4 

be  determined.  The  analyst  is  further  requested  to  specify  the  amount  by 

which  the  conditional  mean  of  x2  will  change  when  x^  deviates  from  ^ 

by  some  fixed  amount.  The  ratio  of  these  two  items  (change  in  ^^yamount 

x  deviates)  will  be  called  m  .  Finally,  the  analyst  will  be  asked  to  state 
1  i! 

the  amount  by  which  the  conditional  standard  deviation  of  x^  will  change 

when  x^  deviates  from  The  ratio  of  these  items  (change  in  g^y/ 

amount  x.  deviates)  will  be  called  m  . 

1  o 

The  resulting  relationships  become: 


K  =  m 
a  H 


E1  =  0 


Kc  ■  A1  -  V. 


E  =  m2 
2  a 


E 

“3  "  3 

*b 


E  =  ~-L 
*b 


where  the  E^  terms  are  the  additive  moments  of  e . 

In  the  special  case  where  m  is  believed  to  be  zero,  the  following 

a 

relation  should  be  used: 


u  +  K 
H2  a 


(xi  ■  ■*)■ 


In  this  case. 


=  A.  K  =  m  E  =  0  E„  =  A. 


E„  =  K 
3  3 


E4  ’  A4  • 
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If  the  analyst  prefers  to  specify  ^the  ratio  of  the  change  in  the 
ional  mode  of  x.  to  the  d< 

mt 

following  equation  may  be  used: 


V 

conditional  mode  of  xg  to  the  deviation  of  x^  from  ^ 


Vo  m  ,  the 

y  u 


rn  =  +  & 

u  x 


m 

_a 


where 

A  is  the  difference  between  the  conditional  mean  and  mode  of 
x2  ,  given  x1  =  ^  . 


More  General  Form  of  Linear  Dependence 


Let  x^  ,  xg  5  .....  x^  be  random  variables  which  are  linearly 

dependent  on  the  independent  random  variables  a^,  a2> .  a^,  It  is 

desired  to  express  the  sum 
n 

l  v 

k=l 

as  the  sum  and  product  of  independent  random  variables.  If  each  x  can  be 
expressed  as 

m 

\  “I  (V^r^y0*:  . 

where 


\i  ■  V  “d  c 

UJ 

%! 


constants; 

the  mean  of  a.  ;  and 
J 

a  random  variable  of  zero  mean  and 
independent  of  all  a' s  and  other  €'s; 


then  the  sum  may  be  expressed  by  the  following: 


which  is  the  desired  result. 

Dependence  Among  Additive-Multiplicative  Combinations  of  Independent 
Elements 

The  situation  sometimes  arises  in  which  independent  random  variables 
reappear  in  more  than  one  expression.  A  simple  example  is  as  follows: 


X  =  X  X 

T  12 


+  X1X3  +  X1X4 


where  the  x^  terms  are  independent  random  variables. 

The  terms  in  the  above  example,  therefore,  are  not  independent. 
This  situation,  however,  may  easily  be  handled  by  factoring  x^  so  that 


x„,  =  x. 
T  1 


x  +  x„  +  x. 
2  3  4 


More  complex  situations,  however,  such  as 


x„,  =  xx  +  xx  +  x,  x„ 
T  1  2  2  3  1  3 


cannot  be  dealt  with  as  simply  since  complete  factorization  ^ 10 1  is  not 
possible.  In  this  case  it  is  necessary  to  take  advantage  of  the  fact  that 
the  expected  value  for  the  sum  of  random  variables,  whether  or  not  they 
are  dependent,  may  be  computed  by  adding  their  expected  values.  The 
reader  will  recall  that  the  k^1  multiplicative  moment  for  x^  is 


XT  =  ^  *T  =  E  4  • 

where  E  designates  the  expected  value. 
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Hence  , 

k 

M.  f x  =  E  /x,x.  +  xx  +  x, x A 

HTJ  V1  2  2  3  1  V  * 

The  preceding  expression  may  be  evaluated  by  expanding  the  argument  and 
then  substituting  for  the  origin  moments  involved.  For  example, 

„  /  2  2  2  2  2  2  2  2 
x_  =  E  { x,  x„  +  xn  x„  +  x,  x„  +  2x,  xx  +  2x  x,x„ 

O  T  VI  O  O  O  1  O  1  2^  V  1^ 

+  2x3\x2') 

=  M2(Xl)M2(X2)+  M2(X2)M2(X3>  M2(Xi)M2(X2) 

+  2M2CXl)“l(X2)Ml(X3) 

+  2M2(X2)  ^Cl) 

+  2M2(xa)  Mi(xi)  mi(x2) 

This  process  is  most  readily  executed  by  automated  means. 
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APPENDIX  II 


PRODUCT  OF  INDEPENDENT  RANDOM  VARIABLES 


If  x  and  x„  are  two  random  variables  with  a  joint  probability  density 
function  of  P^  and  z  represents  some  function*  f  x,(  x^  ,  of 

these  two  variables,  the  probability  density  for  z  is 


m- J  pi2k{f.  X2V2 

“  CD  L 


1  dx„ 


where 


(Z'X2)  = 


=  the  function  that  results  from  solving  f 


x^  in  terms  of  z  and  xg  ;  and 


(xr  x2) 


(^Xl/Sz)  =  derivative  of  %  (^Z»  xy  holdinS  x2  constant. 

In  the  case  under  consideration,  f  x^  =  x^  •  x^,  p12  x^  = 

pi  ©  p2  (<)  •  ei  (*•  x3)  =  z/x2  •  n”ls 

r+co  -i 

P(Z>=  J_w  P1  (z/x^P2(>)X2  ^2  • 

If  =  P2(^2)  5  0  when  their  arguments  are  either  less  than  so 

non-negative  number,  a,  or  greater  than  a  +  1,  the  distribution  for  z 
becomes 


*  This  function  must  be  differentiable  and  single-valued  over  the  ranges  of 
x^  and  x2  where  p12  (xj_,  x£j  is  nonzero.  If  the  inverse  function 
gl  (z,  x^\  is  not  single-valued,  then  the  integral  for  p(z)  must  be 
appropriately  broken  into  segments  and  the  results  added. 
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ru  i 

Pm<’>'  JL  *2  P1  C/XJ>  P2  5)  *2  1 

where 

L  =  a,  u  =  z/a  when  z^  a  (a+1)  ;  and 
L  =  ’  U  =  a+1,  when  z  ^  a  (a+1)  . 

n  0 

and  a  .£  z  £  (a4 1)  . 

If  we  now  define  the  variables  and  functions  as  follows: 

Sty'll  a+1 

In  ^  In  -Tj-  +  in  a,  In  =  Ugln  -jp  +  In  a,  In  z  =  w 
.  a+1 

In —  +  2  In  a 
a 

pi(xD'gi(ui>  p2<fi>"e2<3«>  y,,=g«M  ■ 

it  follows  that 

U 

CW'Ui)  e2  Ca2>  dU2  • 

il 

u 

where 

L  =  0,  U  =  w  when  w  <■  1  ;  and 

u  u  s 

L  =  w-1,  U  =  1  when  w  •>  1  . 

u  u 

The  ranges  of  the  variables  are:  0  ^  u,  ^  1, 0  <  u  ^  1,  and 

x  4 

0  ^  w  ^  2. 

If  u  and  u  are  considered  as  independent  random  variables, 

X  mt 

then  w  =  u.  +  u  .  This  relationship  was  utilized  in  Section  IV  when 
1  3 
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computing  the  exact  solution  for  those  illustrative  examples  that  involve 
products. 

It  becomes  evident  from  thecentral  limit  theorem  and  the -latter 
relationship  that  the  limiting  form  as  n  _  „  of  the  distribution  for 

n 

TT 

k=l 

(where  tne  terms  are  independent  random  variables,  with  a  common 
distribution  defined  for  positive  values)  is  log  normal.  The  log  normal  is, 
of  course,  a  skew  unimodal  distribution. 
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APPENDIX  m 

MOMENTS  FOR  SUMS  AND  PRODUCTS  OF  INDEPENDENT  VARIABLE'S 


The  properties  ascribed  to  the  additive  and  multiplicative  moments 
defined  in  the  glossary  will  be  demonstrated  below* 

ADDITIVE  MOMENTS 

The  additive  moments  A^  ,  A^  ,  ,  and  A^  for  the  sum  of  inde¬ 

pendent  random  variables  x  and  x  are  defined  as  follows: 

1  i 


where 


<k) 

(k) 


m  811(1  n2  = 

E 


the  k^1  origin  moment  of  the  sum 
the  central  moment  of  the  sum 
the  means  of  x  and  x  ,  respectively 

1  it 

the  expected  value 
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Evaluation  of  the  right-hand  members  of  the  above  equalities  yields 


*&s>  -• 

2  2  2 
E (*iVi+x2~vi)  =  E (V^i)  +2E(V^)E(V^)+^V^  ’ 

2  2 

=  E(xrt)  +E(v^)  ■ 


_  E  +  3E(Y“i)  e(Vw5> 

+  3E  E(Ji't‘i)3- 

*  E(v»i).  +  e(VmJ)3  ’ 


E  (YYY^ 4  -3E2(xi-wi+Vtjj)2‘E(^riii)4 

+ 4E  ^i-ki)3  E(v^)  ’ 
+6E(v»i)2e(v^2  ■ 
+4E(xrt)E(v^3 
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It  is  obvious,  therefore,  that  the  multiplicative  moment  for  a 
product  of  independent  random  variables  equals  the  product  of  the  multi¬ 
plicative  moments  of  the  variables, 

CONVERSION  BETWEEN  MOMENT  SETS 

The  following  formulas  show  additive  moments  may  be  generated 
from  multiplicative  moments,  and  vice  versa: 

M1  =  Ar  M2  =  A2  +  A1  ’  M3  =  A3  +  3A1A2  +  A1  ’ 

M4  -  A4  +  3A/  +  4A1A3  +  6Al\  +  Al4  * 

A1  =  Ml*  A2  =  M2  '  Ml2'  A3  =  M3  "  3MiM2  +  2M13 

A  =  M  -  3A  2  -  4MiM  +  6M  2M0  -  3M  4  , 

4  4  2  1  3  1  2  1 
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APPENDIX  IV 

DETERMINATION  OF  THE  PARAMETERS 
FOR  THE  APPROXIMATING  LINEARLY  SCALED  BETA 


Given  the  first  four  additive  moments,  a  linearly  scaled  Beta  distribu¬ 
tion  is  uniquely  defined.  To  determine  the  four  parameters,  a,  g,  a,  and 
b  (see  Glossary),  the  following  procedure  may  be  used, 

1)  Solve  the  following  two  equations  simultaneously  for  a  and  g: 

2 

A3  4(a  -  92  (a  +  g  +  3) 

A2  (a+1)  (g+1).  (a+gH)2 

A4  ,  c  _  3(Q*g+3)  [2  (a+y-2)2  +  pH)  (g+l)  (a+g-4)] 

A2  "  (a+^3)  (g+l)  (cc<-gf4)  (a+g+5) 

U 


where 

a.  <  3  when  A  >0; 

O 

a  -  g  when  A  =  0;  and 

u 

a  >  gwhen  A  <0, 

u. 

2)  Substitute  a  and  3  into  the  equations  shown  below; 
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where 


A 


APPENDIX  V 

ILLUSTRATIVE  EXAMPLES 


INTRODUCTION 

The  examples  shown  in  this  appendix  demonstrate,  under  a  variety  of 
situations,  the  integrity  of  the  approximating  procedures  (see  Section  II) 
built  into  the  computer  program.  The  results  of  the  procedure  are  compared 
with  the  results  obtained  from  actual  performance  of  the  convolution  operation 
for  the  sum  of  independent  random  variables  and  with  the  operation  described 
in  Appendix  n  for  the  product  of  independent  random  variables.  The  latter 
operations  were  performed  on  an  analog  computer. 

EXAMPLE  ONE 

Four  independent  random  variables  describing  the  uncertainty  in  cost 
elements  all  have  the  same  distribution.  This  distribution  can  be  represented 
by  a  linearly  scaled  beta  with  the  following  parameters:  the  most  likely  cost 
is  100  thousand  dollars,  the  most  optimistic  cost  is  50  thousand  dollars,  the 
most  pessimistic  cost  is  250  thousand  dollars,  and  the  central  range  is  110 
thousand  dollars.  The  uncertainty  for  the  total  of  these  four  costs  can  be 
examined,  A  plot  of  the  probability  distribution  for  the  total  cost  is  presented 
in  Figure  7. 

EXAMPLE  TWO 

Two  independent  random  variables,  describing  the  uncertainty  in  cost 
elements,  have  the  same  distribution,  which  can  be  represented  by  a  linearly 
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Figure  7.  Plot  of  the  Sum  of  Four  Independent  Random  Variables 
with  Beta  Parameters 
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scaled  beta.  The  probability  distribution  for  the  sum  of  the  two  costs  is 
desired. 

Four  sets  of  parameters  for  the  linearly  scaled  beta  are  considered. 
(The  values  for  each  parameter  are  shown  in  thousands  of  dollars.) 


Optimistic  (X.^)  ^reliable  (Xp)  Pessimistic  (X^)  ^lage**fe  ) 


50 

100* 

150 

64 

50 

100** 

150 

19 

50 

100*** 

250 

110 

50 

100**** 

250 

38 

The  distributions  for  total  cost  appear  in  Figures  8  to  11. 

EXAMPLE  THREE 

Two  independent  random  variables,  describing  the  uncertainty  in 
planning  factors,  have  the  same  distribution  which  can  be  represented  by  a 
linearly  scaled  beta.  The  most  pessimistic  value  for  each  of  the  random 
variables  is  200  percent  of  the  most  optimistic  value.  The  probability 
distribution  for  the  product  of  the  two  costs  is  desired. 


*The  normalized  form  of  this  distribution  appears  as  Curve  III  in 
Figure  1. 

**The  normalized  form  of  this  distribution  appears  as  Curve  I  in 
Figure  1. 

***The  normalized  form  of  this  distribution  appears  as  Curve  m  in 
Figure  2, 

****The  normalized  form  of  this  distribution  appears  as  Curve  I  in 
Figure  2. 


100  200  300 

THOUSANDS  OF  DOLLARS 


Figure  8.  Plot  of  the  Sum  of  Two  Independent  Random  Variables 
with  Beta  Parameters 
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Figure  9.  Plot  of  the  Sum  of  Two  Independent  Random  Variables 
with  Beta  Parameters 
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Figure  10.  Plot  of  the  Sum  of  Two  Independent  Random  Variables 
with  Beta  Parameters 


Figure  1  ] .  Plot  of  the  Sum  of  Two  Independent  Random  Variables 
with  Beta  Parameters 


Four  sets  of:parameters  for  the  linearly-  scaled  beta  are  considered. 


Most. 

Central 

Optimistic  OH 

-  ,±j 

Probable  (Xp) 

Pessimistic  iX^y 

Range  (( 

'  1.  ' 

1.5* 

2- 

.64 

1 

1,5** 

2 

.19 

1 

1.-25*** 

2 

.55 

i  . 

1.25**** 

2' 

.19 

The  distributions  for  the  products  appear  in  Figures  12  to  15. 

EXAMPLE  FOUR 

This  example  is  similar  to  Example  Three  with  the  exception  that  the 
most  pessimistic  value  in  this  case  is  120  percent  of  the  most  optimistic 
figure.  Four  sets  of  parameters  for  the  linearly  scaled  beta  are  considered: 


Optimistic  (X^) 

Must 

Probable  (X^) 

Pessimistic  fX^l 

Central 
Range  (( 

5 

5.5* 

6 

.64 

5 

5.5** 

6 

.19 

5 

5.25*** 

6 

.55 

5 

5.25**** 

6 

.19 

The  distributions  for 

the  products  appear  in  Figures  16  to  19, 

♦Normalized  form, 
♦♦Normalized  form, 
♦♦♦Normalized  form, 
♦♦♦♦Normalized  form, 


Curve  HI  in  Figure  1. 
Curve  I  in  Figure  1, 
Curve  HI  in  Figure  2, 
Curve  I  in  Figure  2. 
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Figure  12.  Plot  of  the  Product  of  Two  Independent  Random  Variables 
with  Beta  Parameters 


Figure  13.  Plot  of  the  Product  of  Two  Independent  Random  Variables 
with  Beta  Parameters 
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1.00  1.15  1.30  1.45  1.60  1.75  1.90  2.05  2.20  2.35  2.50 

Figure  15.  Plot  of  the  Product  of  Two  Independent  Random  Variables 
with  Beta  Parameters 


Figure  i6.  Plot  of  the  Product  of  Two  Independent  Random  Variables 
with  Beta  Parameters 
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Figure  18.  Plot  of  the  Product  of  Two  Independent  Random  Variables 
with  Beta  Parameters 
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25.0  25.5  26.0  26.5  27.0  27.5  28.0  26.5  29.0  29.5  30.0 


Figure  19.  Plot  of  the  Product  of  Two  Independent  Random  Variables 
with  Beta  Parameters 
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APPENDIX  VI 


DESCRIPTION  OF  THREE  7090  FORTRAN  PROGRAMS 

INTRODUCTION 

This  appendix  describes  a  series  of  three  7090  FORTRAN  programs 
(see  Figure  20)  which  convert  a  cost  analyst's  subjective  feelings  about 
uncertainty  in  cost  estimates  for  system  elements  into  quantitative  data  to 
express  Hie  uncertainty  in  the  total  system  cost, 

DESCRIPTION  OF  THE  PROGRAMS 
Program  I 

This  program  (see  Figure  21)  accepts  as  inputs  the  four  parameters 
specified  by  the  analyst  for  each  element,  and  from  these  computes  a  linearly 
-scaled  beta  distribution  which  approximates  the  distribution  described  by  the 
parameters  for  die  cost  of  the  element.  This  distribution  is  presented  in 
terms  of  the  parameters  describing  the  shape  of  the  curve  and  the  first  four 
additive  moments  of  the  distribution. 

The  four  input  parameters  are: 

(a)  XL  -  die  lowest  cost  (lower  1  percent  tail); 

(b)  Xg  -  the  highest  cost  (upper  1  percent  tail); 

(c)  Xp  -  the  most  probable  cost  (mode);  and 

(d)  CR  -  the  length  of  the  80  percent  central  range 

(upper  10  percent  tail  minus  lower  10  percent  tail). 


|  PUNCH  OUTPUT 
SET  UP  CAROS 
(NSERTfiN 
PROGRAM  3C 


PROGRAM  31 

r — T  “ 


PERFORM 

OPERATIONS 

ON  MOMENTS 


PRINT  OUT 
OUTPUT 


INSERT  PUNCHEO 
CAROS  OUTPUT 
INTO  PROGRAM  3ZC 


PROGRAM m 


Figure  20.  Block  Diagram  for  Programs  1,  II  and  1 11 
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•  .*v- 


READ  NEXT  INPUT  CARD 
COMMENT  XL,  XH.XP.  CR 


r END  OF  DATA 
CARDS  CALL 
EXIT 


>190 


YES 


cr>o; 


WRITE  OUT 
COMMENT,  SET  UP 
FORMATS  ft  WRITE 
OUT  INPUT  DATA 


•178'.., 


■JLL. 


COMPUTE. 

X  .JSt— 

*  XH-XL 
CRN  «CR/XH-XL 


NEGATIVE 

660-263 


300-182 


Of.fi  K  *  1 

A  , 

8<r2 

,«ar+i 


fH.u» 


-  3/2 


VTU) 


«:»)+/» (11+2; 
CO  T0:;i6:;  . 


POSITIVE; 

MAX  CR  EXCEEDED 

GO  TO  33 

COMPUTE: 

A  *  0.H8-603«u  +  3.96  »u2 
ALF  «  L76- 3;9* (1+29,16*1)* 

<r  .  0:39  •  CRN/e  L?A»CRN**'AL  F 
K«  XP-XL-/XH-XP 
SET  UP  NEEDED  CONSTANTS! 
V^..  (QlHfl-H) 


POSITIVE 


CALL  CU8ERT 


4-172 


HREALS^IS  ONLY 
I  REAL  SOLUTION 
fill)  •QREAL3 

<a/<l)Vj8(l)*R 

*IND.2*1. 


WE  NEXT  TEST  TO  S' 
I  SEE  IF'ONE  REAL  ' 
’  SOLUTION  OR 


ITHREE  IS  RETURNED! 


3- 166 

QREAL1  ISA 
REAL  SOLUTION 
<fi  U)i  OREAU 
«  TU  *  fiWik 
TN02--L  . 


7-177 


|;0REAL2  IS  ONLY 
.  REAL ROOT 
■fi  (1)  *0RCAL2 
•  (U*/9tl)*R 
IND  2  «1 


THERE  ARE  3 'REAL 
ROOTS  SET 
£(2):>  0RE  AL2 
a  (2)»0<2)*k 
fi  (3)  ■  0REAL3 
O  13)  i  filiUk 
IND2  •  3  . 


{COMPUTE ) 

Figure2i.  Program  Is  MainRoutine 
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.4(«-0H«+/3  +  3)  . 

(a+lM/S+l)  (a+#  +  4  ) 

9(a+/9+3)£t  (•+£  +  -2  )*  +  {a+l)(0+i){a+0-4)] 
(o  +  U(S+l>(o+/5  M)(a>0  +5) 


NOTE:  The  cential  range  may  not  exceed 

>•&**£)  - °.«  I  <«(yxi)  -5»-\)  I- 

If  it  does,  the  program  will  skip  the  case  and  go  on  to  the  next  case.  The 
beta  distribution  is  computed  as  follows.  First,  the  standard  deviation  a 
is  determined.  Compute 


~  .  VXL 

'  VXL  ' 

1 

c  -  °B 

\  W 

Let 

U  =  |  0.5  -  X  j 

Then 

0.39  0, 

.  rn 

C  -  a 

1  +  ACtt 

rn 

where 

A  =  0.148  - 

0.603U  +  3.98U2  ,  and 

a  =  1.78  -  3.9U  +  29.16U2  . 

These  equations  were  determined  by  curve-fitting,  using  several  different 
values  for  the  variables  involved.  Then  determine 


v* 

vs- 
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The  parameters  a  and  p  are  then  determined  by  finding  the  real  solution 
of  the  following*  pair  of  equations  (see  Appendix  II): 


K  =  c/P  ; 

g2  _  (q+i)  rn) 

(a-p+2)2  (ctfP+3) 


For  each  pair  of  roots,  compute 


n  = 


a+i 


Let 


a+p+  2  * 

~  VXL 


X  = 


a  . 


xH-xt  m 


and  find  all  a,  $  which  satisfy 

a  -  0.5 

o  <  c..— Yas.  <  i 

W  M  <V  _  *•  • 

X-  0,5 

(In  general,  there  Will  be  only  one  such  pair.) 

When  this  pair  (a,Q  is  determined,  compute : 


4  (cup)  (ctf$f3) 

“i  =  T  * 

(a+l)  (PH)  (afpf4) 


*The  equation  is  solved  in  the  equivalent  form: 
o2  (K+l)3  9  +  [7o2  (K+l)2  -  k]  p2  (K+l)  p  +  1202  -  1  =  0  . 
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and 


*2  = 


3  (ch#3)  [  2  (ctf$f2)2  +  (otfl)  ($t-l)  'otWJ-4)] 
(a+l)  (fri)  (oa-M)  (»0i-5) 


All  necessary  parameters  are  thus  computed  for  determining  the  four 
moments.*  These  are  computed  as  follows: 

(a)  first  additive  moment  =  mean  =  X  = 

(b)  second  additive  moment  =  (standard  deviation)2 = Y 2 

(c)  third  additive  moment  =  V  -  +-^  Y2  . . . .  ** 


]^<J 


■>-.sign 


;  and 


(d)  fourth  additive  moment  =  W  =  ^  ~^Y2  * 

Also  included  in  the  output  are  a  g  and  ~\f?  (the  standard  deviation). 


Program  fl 

The  analyst  now  has  cost  distributions  (summarized  by  four  moments) 
for  each  of  the  elements  of  the  system  under  consideration.  He  specifies  the 
additive-multiplicative  combinations  of  these  distributions  that  Will  yield  the 
total  system  cost,  and  inserts  the  moments  of  each  distribution  into  Program 
H  (see  Figure  22).  Program  H  computes  the  specified  combinations  and 
prints  out  the  moments  of  the  total  cost  distribution,  as  well  as  certain  other 
subtotals  of  interest  to  the  analyst. 


*In  the  main  text,  the  four  additive  moments  are  referred  to  as  A^^  ,  A  , 
A3  ,  and  A4;  in  this  appendix  they  are  referred  to  as  X,  Y,  V,  and  W, 
respectively. 
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Figure  22.  Program  11:  Main  Routine 
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Additive- Multiplicative  Combinations  of  Independent  Random  Variables 
Independent  random  variables  D.  and  D  are  given;  the  distributions 

X  u 

for  each  summarized  by  its  first  four  moments.  Then 

D1  =  D1  (\»  Yx»  Vx»  W^)  »  “d 

D  =  D„  (X.  Y  ,  V  ,  W  'N  . 

2  2\^  2  2  2  2/ 

The  distribution  for  D  +  D  is  computed  by  component  adding 

X  it 

moments  on  the  basis  of  components;  that  is, 

VVVVV1^  • 

The  distribution  for  D  D  is  computed  by  the  following  equations:* 

X  2 

=  Dm  (X,  Y,  V,  W) , 

where 

X  =  X,Xn  . 

1  2  * 


Y=YiVxi2 

Vx22V 

)+VlX23+V2X13  +  6XlYlX2Y2i“‘d 

ViW3C 

W  =  WXW2  +  3  ( 

"W1Y22+W2Y12) 

+  4(W1X2V2  +W2X1Y0 

+  6 

+’Y^Y^\  +  WX^+WX4 

1  1  1  2  J  W12  W2  1 

♦This  is  the  result  of  solving  simultaneously  the  equations  shown  in 
Appendix  III. 
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These  formulas  are  equivalent  to  converting  tine  additive  moments  into 
origin  moments,  multiplying  these  origin  moments  in  terms  of  components 
and  reconverting  the  product  moments  to  additive  moments. 


Additive-multiplicative  combinations  of  distributions  may  be  computed 
by  evaluating  a  composite  term  as  in  an  algebraic  expression.  For  example, 
to  evaluate  ^  +  D^D ^  ,  first  evaluate  DgD4  ,  add  D 2  to  the  result, 
and  multiply  the  last  result  by  .  NOTE:  This  expression  may  not  be 
evaluated  in  the  form 


D  D  +  D  D  D 
12  13  4 


since  the  two  product  terms  are  not  independent  (both  depend  On  D^)  . 


Additive-multiplicative  combinations  of  distributions  are  specified  for 
the  program  on  two  levels:  lines  and  groups  (see  Figure  23). 


lines 

A  line  is  defined  as  an  additive- multiplicative  combination  of  distribu¬ 
tion  of  the  following  form: 

AAA  {AA  [AA  (0+A0)  +  AA0  +  AO  +  o]  +  AO  +  AO  +  AO  +  0} 
123  4  5  6  7  8  9  10  11  1213  14  15  16  1718  19  20  21  22  23 

where  each  "A"  or  "0"  symbol  represents  a  distribution  specified  by  its 
first  four  additive  moments  X,  Y,  V,  W,  The  distribution  for  "A"  is 


LINES 

GROUP  TOTALS 

LINES 

GROUP  TOTALS 

LINES 

GROUP  TOTALS 

Figure  23. 


Generalized  Cost  Structure  Built  into  Program  )] 
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defined  by  the  moments  (1,  0,  0,  0),*  for  "0"  by  the  moments  (0,  0,  0,  0) 
It  may  be  observed  that: 


(a)  if  "0"  is  added  to  any  distribution  ,  the  resulting  moments 

are  simply  those  of  the  distribution  ;  that  is,  0  +  ; 

(b)  if  "0"  is  multiplied  by  any  distribution  ,  the  resulting 
moments  are  (0,  0,  0,  0);  that  is,  =  0;  and 

(c)  if  "A"  is  multiplied  by  any  distribution  ,  the  resulting 

distribution  moments  are  simply  those  of  that  is,  AD^  =  . 


In  other  words,  "0"  and  "A"  resemble  the  0  and  1  of  ordinary  arith¬ 
metic,  Because  of  this,  the  line  total  operator  as  it  stands  is  equivalent  to 
0,  If  any  "0"  were  replaced  by  some  D^,  the  line  total  would  become  D^. 


Lines  are  specified  to  be  equivalent  to  certain  additive- multiplicative 
combinations  of  distributions  by  the  following  technique:  the  analyst  is 
permitted  to  replace  any  term  in  the  line  total  operator  by  whatever  distribu¬ 
tion  he  desires.  Such  substitution  will  alter  the  value  of  the  operator.  For 

th  tii 

example,  by  replacing  the  7  ,  8  and  10  terms  with  the  distribution  , 
D2  ,  and  D3  ,  respectively,  the  line  will  be  made  equivalent  to  the  combi¬ 
nation  (p2  +  D2^),  This  can  readily  be  seen  by  representing  the  operator 
symbolically: 


L11  jl*l  l'D^D^l'D^  +  1*  1-  0  +  1-  0  +  0  +  1*0  +  1'0  +  1*0  +0j>  , 


Application  of  the  obvious  rules  reduces  to  this 


1 


[Dl(D2  + 


d3)+(»]^  =  d1(d2  +  d3) 


♦The  slashed  zero  (0)  is  used  to  denote  the  numerical-zero  to  avoid 
confusior  , 
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Therefore,  simply  by  specifying  the  iine  and  the  term  into  which  a  given  set 
of  distribution  moments  are  to  be  substituted,  a  large  variety  of  combina¬ 
tions  can  be  evaluated. 

Groups 

On  the  second  allowable  level  of  combination  of  distributions,  the 
various  line  totals  are  separated  into  groups,  corresponding  to  subsystems 
of  toe  system  under  consideration.  First,  all  line  moments  for  a  group  are 
Computed,  Then  these  moments  are  summed  by  components  to  produce  a 
term  S.  This  term  may  be  combined  with  other  distributions  in  toe  following 
operators: 


(Gl) 


AS(A  +  AO) 
1  2  34 


,  and 


(G2)  AS(AO) 
1  34 


Combinations  are  specified  in  toe  same  manner  used  for  lines:  any  A  or 
0  may  be  replaced  by  toe  moments  of  any  distribution;  A  term  replaced  in 
toe  (Gl)  operator  is  also  replaced  in  toe  .(G2)  operator  by  toe  same  moments 
except  for  the  second  term,  which  appears  only  in-toe  (Gl)  operator,  it  may 
be  noted  that  if  no  term  is  replaced  in  toe  operators;  toe  result  of  (Gl)  is  S 
and  toe  result  of  (G2)  is  "0,” 

Final  Total 

All  moments  obtained  from  toe  (Gl)  operator are  added-by  components 
to  give  toe  final  total.  Thus,  toe  (Gl)  group  total  shoUid  be  sO  defined  that 
the  total  system  cost  distribution  is  toe  sum^of  all  Of  these  totals. 
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Program  in 

From.  Program  II,  the  analyst  has  moment  sets  summarizing  distribu¬ 
tions  for  subsystems  and  total  system  cost  uncertainty.  These  moment  sets 
are  now  inserted  into  Program  m  (see  Figure  24),  which  converts  the 
moments  into  the  beta  distribution  parameters  and  the  bounds  for  the  curve 
and  for  any  specified  percentage  tails. 


From  the  input  moments,  X,  Y,  V  and  W  for  each  case,  determine 


The  followinjg  equations  are  now  to  be  solved  for  and  g  ,  the:beta 
curve  parameters: 


1  (on).on)  (»m)2 


j.  .  3  il  2  2  foe- ii  mv 

2  "  .  (cu-1)  (gflj  '(atm)  -(wff-5) 


The  constants  g^  and  ^  determine  whether  or  not  there  is  a  mean¬ 
ingful  solution'  to  these  equations.  *  (Refer  to  Appendix  H  for  the  details  of 
toe  solution  to  these  equations. )  if  there  is  a  solution  to  these  equations, 
use  a  g  to  compute  for  toe  normalized  beta  curve: 


a+1 

HN  =  ~~r  ■  (normalized  mean)  ; 

crp* 


*For  a  unimodal  beta  distribution,  a  and  g  must  be  positive  and  real. 
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»?- 


■'  •  r&G3Z 


Afern  <»i) 

N  (Ort-^2)  ~\J(»&r3 


(normalized  standard  deviation);  and 


X 


a 

a+p 


(normalized  mode)  . 


The  normalized  curve  must  be  integrated  to  determine  the  total  area 
under  the  curve:  Integrate 


where 

bL  =  UJJ  -  4  ! 

bH  ’  Bto(1*  MN  +  4V  • 


This  integration  is  performed  numerically  with  Simpson's  rule.  The  number 
of  intervals  used  is  NINT,*  an  input  (800  is  suggested).  Determine  the 
desired  area  under  each  tail: 


A  A  =  A  x  AREA, 


where 

A  is  an  input  specifying  the  fractional  tail  desired. 


*Number  of  intervals  to  be  used  in  the  Simpson  integration. 
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The  curve  is  integrated  again,  one  interval  at  a  time,  from  each  end, 

* 

until  the  X-values  which  are  the  bounds  for  the  tails  are  reached.  These 
are  called  XNCL  and  XNCH,  respectively.  From  the  bounds  for  the  normal¬ 
ized  curve,  the  bounds  for  the  scaled  beta  curve  are  then  determined. 
Compute: 


jr-X 

(mehn); 

a  -  ~\Jy 

(standard  deviation)  . 

Then 

the  lower  bound  is 

aPN. 

XL 

the  upper  bound  is 

301  ■  • 

the  lower  tail  bound  is 

XGL  ^XNCL  -  ^  ;  and 

the  upper  tail  bound  is 

XGB  (X.NCH-.,n)  . 

USE  OF  Tilt  PROGRAMS 
Program  1 


The  input  for  Program  I  consists  of  the  four  parameters  specifying 
the  uncertainty  in  the  cost  of  each  individual  system  element: 

X_  =  the  lowest  cost  {lower  1  percent  tail); 

L» 

Xj.  =  the  hipest  cost  (upper  1. percent  tail); 
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i 


the  most  probable  cost  (mode);  and 
the  length  of  the  80  percent  central  range  (difference 
between  upper  and  lower  10  percent  tails). 

The  data  for  each  element’s  cost  uncertainty  is  punched  on  a  separate 
input  card.  Each  card  has  the  following  format  (refer  to  Figure  25): 

■Column  1  is  blank. 

Columns  2  to  24 contain  a  Hollerith  -identification  which  the  user  may 
specify  to  identify  the  data.  This  identification  will  be  printed  out 
before  the  data  for  the  case  specified  by  the  card  when  the  program 
output  is  presented. 

,  Columns  25  to  32  are  left  blank. 

and  CR. 

The  data  must  be  in  E-format,  containeddn  a  field  of  ten  columns. 
The  first  column  of  each  field  should  be  blank  (columns  33,  43,  53, 
63).  The  second  contains  the  sign  (if  positive,  it  may  be  omitted). 
The  third  through  sixth  columns  contain  the  three  most  significant 
digits  of  the  number  in  the  form  X*XX,  The  seventh  through  tenth 
columns  contain  the  power  of  10  by  which  the  significant  digits  must 
be  multiplied  to  give  the  actual  number.  The  exponent  must  be  in 
the  form  E±XX  (+  sign  may  be  left  as  blank).  Thus,  the  entire 
number  is  shown  by: 

±  (X)  .  (X)  (X)  E  ±  PQ  (^ 

IMPORTANT:  The  last  significant  digit  of  the  exponent  must  be  in  the 
tenth  column. 


Columns  33  to  72  contain  the  four  parameters  X^,  X^,  Xp, 


V 


rT' 
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The  cards  may  be  arranged  in  any  desired  order;  the  cases  will  be 
processed  in  the  order  of  their  appearance  in  the  deck.  Following  the  card 
for  the  last  (physical)  case  a  special  card  with  -1.  j#E  ±  ??  punched  in 
columns  44  to  53  must  appear.  This  signals  that  all  data  has  been  processed, 
and  the  program  will  exit. 

Program  Deck 

■"V*  program  deck  should  have  the  following  sequence  of  cards: 

m  job 
in  FMS 

*  XEQ 

[  BINARY  PROGRAM  DECK  ] 

*  DATA 

DATA  CARDS  j 
CR  =  -1.0  card  j 
m  ENfiffOB 

Output 

The  output  data  (see  Figure  26)  will  be  arranged  in  the  following 
manner.  The  identification  will  be  printed  at  the  head  of  the  data  for 
e'.ch  case.  The  input  data  will  then  be  printed;  following  this  will  come  the 
output  data: 

X  =  the  first  additive  moment  (the  mean)  of  the  approximating  beta 
distribution; 

Y  =  the  second  additive  moment  (variance) ; 

V  =  the  third  additive  moment; 

W  =the  fourth  additive  moment; 
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CASE-EXAMPLE  I  PROGRAM  I 


■***■- •*.•< 


o 

o 

o 

o 

n 

o 

o 

o 

o 

u 

-1- 

<UJ 

<  UJ 

<  UJ 

<  UJ 

«u 

h® 

HS 

h  IO 

h<o 

ho 

~  N 
ffl  ♦ 

CO  • 

UJ  ,1 

m® 

S  . 

CD.  * 

CM 

CM 

N 

O 

O 

5 

o 

-O 

O 

<  1 

XUJ 

<  i 

I  w 

<  l 

XU 

CLIO 

o_  *2 

q-2 

0.5 

CLtf) 

JO) 

-J  io 

JO 

Jo 

<{0 

<  r0 

<CM 

<■— 

*K> 

>2 
£C  K) 

d* 

w- 


o 

in 

K  O 
o.w 


*8 

xo 


o 

o 

Ui 

I  o 
X  o 
K) 


.2 

x§ 


CM 

O 

ill 

*8 


(0 

o 

viii 

a:N 
dm 
s  cn 
o 

hi 

uo  *■ 

—  LI 

*S 

<0  I 

o 


in 

o 

t  m 

o-: 
m  w  - 
o 


u 

k2 

uP 


in 

o 


o 

CM 

UJ- 

® 

.in 

’in 


<r 

O: 

LI 

„;0 
x  o 

N 


w 

o 


«• 

>° 
l_:  LI 
*8- 
s- 


in 

:« 

cxi 

i 


CM- 

to. 

d- 

CM  CO 

o 

Li  CD 

crO  o 

oPf 


LI 


O 

Hi  “ 
?w. 
s-  T 


til 

Li 

u 

u 

CM 

0.  o 

o 

0.0 

m 

to 

CM 

0.0 

IO 

O 

x  o 

CM 

XO 

x  tf> 

xo 

o 

Id 

in 

til 

K) 

Ui 

- 

UJ 

10 

UJ 

<0 

CM 

CM 

CO 

CO 

>CM 

>* 

>® 

>9 

>® 

to 

N 

CM 

1 

V 

‘  CM 

10 

CO 

CM 

CM 

o 

o 

o 

o 

o 

UJ 

IO 

K) 

LI 

IO 

o 

u 

zo 

oa 

0 

u 

IO 

** 

0 

• 

XO 

XO 

xn 

xo 

Id 

CD 

= 

CM 

UJ 

o 

CM 

u 

CO 

5 

CM 

Id 

N 

z 

rO 

a> 

id 

o 

>o 

CM 

CM 

CM 

>•9 

CM 

W® 

CM 

>-3 

*- 

CD 

(0 

in 

CM 

IO 

CM 

< 

O 

O 

O 

O 

o 

O 

Z 

\ 

UJ 

-JO 

CO 

O 

a> 

N. 

Li 

JO 

,tn 

K 

5 

Id 

JO 

IO 

CO 

UJ 

CM 

o 

X  O 

o 

u 

XO 

o 

3 

X  — 

o 

-J 

X® 

o 

UJ 

s 

Li 

CM 

UJ 

X 

CM 

u 

m 

£ 

I 

Id 

£ 

* 

Id 

N 

X  CM 

£ 

x5 

£ 

x« 

? 

xco 

5 

. 

x*r 

“ 

* 

N* 

5. 

<r 

-? 

cp 

o 

a 

Q 

to 

o 

CO 

CO 

(0 

CO 

t 

0) 

H 

K 

<0 

K 

CO 

h* 

o 

CL 

1 

u 

3 

a. 

1 

u 

>- 

2 

1 

u 

H 

3 

3 

a. 

1 

u 

K 

0 

3 

a. 

K 

to 

CL 

i- 

CO 

£ 

i- 

CO 

A 

h* 

CO 

(L. 

H 

3 

O 

< 

o 

z 

•D 

O 

< 

o 

z 

3 

O 

< 

o 

lb 

z 

3 

O 

< 

o 

z 

3 

o 

-77- 


Figure  26.  Scruple  Output"  for  Program  1 


^  !  the  parameters  which  determine  fee  shape  of  the  beta 
g  j  curve;  and 

V7  -  the  standard  deviation 

All  data  will  be  printed  in  E-format  to  three  significant  decimal  digits. 
Program  II 

It  is  important  that  the  user  read  the  previous  discussion  of  Program 
n  for  an  understanding  of  how  to  use  this  program. 

Input 

Each  group  is  processed  separately  by  fee  program,  so  all  cards 
relating  to  fee  same  group  must  be  placed  together.  For  each  group,  the 
following  input  cards  must  be  present  (see  Figure  27): 

The  first  card  in  each  set  of  group  cards  must  have  a  positive  number 
JJ  (preferably  fee  number  of  fee  group)  in  Column  5.  (If  it  is  a  two- 
digit  number,  use  columns  4  and  5.) 

The  next  set  of  cards  specifies  fee  line  data  for  this  group.  There  is 
one  card  for  each  set  of  moments  which  are  to  be  used  in  this  group. 
Each  card  has  an  identification  field  in  columns  2  to  24  as  in 
Program  I, 

Column  25  is  blank. 

Columns  26  and  27  contain  fee  group  number,  J,  for  fee  present 
group.  All  cards  for  one  group  must  be  placed  togefeer. 

Columns  28  and  29  contain  fee  line  number,  I,  in  which  this  particular 
set  of  moments  is  to  be  combined.  The  I's  in  each  group  must  be 
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Figure  27.  (concluded) 


arranged  from  1  consecutively  up  to  IMAX,  which  cannot  be  greater 
than  50. 

Column  30  is  blank. 

Columns  31  and  32  contain  the  term  (column)  number  K  for  which 
the  moments  on  this  card  are  to  be  substituted  in  the  line  total 
specified  by  I.  The  K  term  may  be  any  integer  between  land  23, 
inclusive. 

Columns  33  to  72  contain  the  moments  X,  Y,  V,  W  which  are  to  be 
substituted  for  term  K  in  Line  I  in  group  J.  These  columns  have 
the  same  format  as  on  the  Program  I  data  cards, 

•  The  group  total  cards  follow  the  line  cards.  These  have  the' same 
format  as  the  line  cards,  except  that: 

I  is  always  0  for  a  group  total  input;  and 
K  is  the  term  number  for  which  the  moments  on  the  card  will  be 
substituted  (K  ranges  from  1  to  4,  inclusive)* 

A  card  which  is  blank,  except  for  a  "-1”  in  columns  28  and  29,  follows  the 
last  group  total  card  (or  last  line  total  card  if  there  is  no  group  total  data)* 
This  signals  the  end  of  the  group. 

The  cards  for  all  subsequent  groups  are  arranged  in  the  same  way. 

A  card  which  is  blank  except  for  a.  0  punched  in  column  5  appears  after  the 
I  =  -1  card  for  the  last  group.  This  signals  the  end  of  the  data. 

Program  Deck 

m  JOB 
in  FMS 
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* 


♦ 


ZEQ 

[binary  deck] 

DATA 


JJ  =  1  ]  card 
Line  Total  Data  I  >  0 
Group  Total  Data  I  = 
I  =  -1  card 

w  a 

:  JJ  =  2  card 


m  ENDJOB 


Output 

The  line  total  input  data  for  the  first  physical  group  in  the  deck  appears 
on  the  first  page  followed  by  the  group  totalirput,  if  any,  for  this  group  on 
the  next-page  (see  Figure  28).  Each  entire  card  will  be  printed  out  except 
for  the  JJ  and  I  =  -1  cards. 


The  line  total  moments  X,  Y,  V,  W  for  this  group  are  shown  on  the 
next  page.  In  addition,  three  characters  will  be  printed; 

J  =  the  group  number; 

I  =  the  line  number;  and 
L  =  0,  to  indicate  a  line  total. 


The  X,Y,  V,W  are  interpreted  in  the  same  way  as  for  Program  I,  except 
that  now  they  are  moments  representing  the  distributions  of  random  vari¬ 
ables  which  are  aggregations  of  other  system  element  costs. 

The  group  total  output  for  the  group  appears  on  the  next  page. 
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J  =  group  number 
I  =  0,  to  indicate  group  total; 

L  =  1,  for  the  (Gl)  total;  and 
L  =  2,  for  the  (G2)  total 

All  groups  will  have  their  output  in  this  form.  After  the  last  group 
total  output  the  final  total  will  appear.  Its  characters  are 

J  =  0;  I  =  0;  L  =  1. 

All  output  data  will  be  punched  on  cards  for  direct  insertion  into 
Program  HI. 

« 

Program  HI 
Input 

The  first  card  specifies  the  two  parameters  A  and  NINT  for  the 
integration  process  (see  Figure  29): 

Columns  1  to  5  are  blank. 

Columns  6  to  14  contain  A,  the  fraction  of  the  total  area  desired 
in  each  tail. 

Columns  15  to  35  are  blank. 

Columns  36  to  38  contain  NINT. 

Rest  of  card  blank. 

The  next  cards  contain  the  data  to  be  processed  by  the  program. 
Columns  2  to  24  contain  identification. 

Column  25  is  blank. 
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Figure  29 .  input  Format  for  Program  1 1 1 


Columns  26  and  27,  28  and  29,  30  and  31,  and  J,  I,  and  L  indicate 
the  control  characters. printed  with  Program  H  output.  These  are 
merely  for  the  purpose  of  identification;  they  are  not  used  by  the 
program. 

Column  32  is  blank. 

Cblumns  33  to  72  sX,  Y,  V,  W  are  in  same  format  as  Programs  I 
ahd  n. 

A  card  with  nine  columns  (44  to  52)  containing  a  -l.<WE+ffjf  mu3t 
appear  after  these  cards.  This  signals  the  end  of  the  data. 

Program  Deck 

in  JOB 
m  FMS 

*  XEQ 

[programdeck] 

*  DATA  - 

datacards] 

Y  =  l.ff  card] 

in  endjob 

Error  Returns 

A  number  of  conditions  will  cause  the  program  to  indicate  that  no 
solution  can  be  found  for  the  present  case.  When  this  happens,  the  program 
will  proceed  with  the  next  data,  because  the  distribution  specified  by  the 
moments  for  this  case  cannot  be  converted  by  the  program. 
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Output 

The  output  data  (see  Figure  30)  for  each  case  will  be  preceded  by  the 
identification  field  and  the  input  data.  The  output  will  consist  of: 


♦Refer  to  the  definition  of  the  linearly  scaled  beta  distribution  in  the 
Glossary;  a  =  and  a+b  = 
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appendix  vrt 

DETAILS  OF  FORTRAN  SUBROUTINES 


PROGRAM  I 


Program  J  requires  the  soiution  of  a  cubic  equatiou.  The  routine 
CtIBERT  (see  Figure  31)  is  used  (or  this  purpose. 
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PROGRAM  in 


Subroutine  SETUP 

The  subroutine  SETUP  (see  Figure  32)  determines  whether  a  solution 
exists  to  the  equations  in  Program  III  for  (a ,  g)  and  g^  (<x ,  It 
checks  certain  criteria  on  g_  and  n  to  determine  the  nature  of  the  equa- 

i  it 

t'ons  for  a  and  6  as  shown  below. 

1)  The  program  will  skip  any  case  for  which  is  greater  than  4  or 
g  is  greater  than  9. 

A 

2)  If  &  is  less  than  0.  01 ,  no  iteration  is  necessary.  In  this  case , 
the  following  checks  are  made: 

(a)  If  p  is  less  than  2. 9994,  then  the  final  solution  is 

it 

5f32-9 

a  =  e=2(3Hy  * 

(b)  If  3  is  between  2. 9994  and  3.0015,  the  final  solution  is 

it 

approximately 

a  =  g  =  00. 

These  are  the  answers  used  in  this  case. 

(c)  If  p  is  greater  than  3.0015,  there  is  no  solution,  and  this 

it 

case  will  be  skipped. 

3)  If  g  >0. 01,  g  is  compared  to  3. 01  with  the  following  results: 

X  it 

(a)  If  g  is  less  than  3.  (/(ft,  iteration  is  necessary,  and  procedure  (4) 

it 

is  used. 
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ENTER 


Figure  32.  Program  111;  Subroutine  SETUP 
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(b)  If  is  greater  than  3. 001,  then  the  trial  solution  for  g  is  +  ®. 
In  this  case,  equations  are  derived  for  g^  and  g^  with  g  =  ®, 
and  the  value  of  a  determined  in  each  equation.  The  equations 
are: 

«.<K)  ■  V 1 

9  -  g„ 

a»  (?t)  =  gg  -  3 

If  a  ,  there  is  no  solution  possible,  and  the 

case  will  be  skipped. 

If  a,  ,  then  the  number  J  gg  g^^  ,  g=“J 

-  g  I  =  beta  is  computed.  Should  g  <  0. 0025  and 

Ct  | 

4  8 

a  10  »  then  this  a  g,  is  the  solution  for  a,  and  10 

oo  Vly-  co  j. 

is  used  for  g. 

If  either  of  these  conditions  is  not  satisfied,  iteration  is  used, 
and  the  program'proceeds  according-to  the  next  step. 

4)  When  iteration  is  used,  the  initial  guess  for  a  and  g  is  deter¬ 
mined  thus: 

(a)  If  g  <  2. 929,  then  the  guess  for  g  is 

5g2-9 

S  =  2lp2)  ‘ 

a  =  g/2  is  used  for  the  RAFSON  routine  and  a  =  g  for  REDOIT. 
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(b)  If  80  >  2. 929,  then  a  =  8  =  50  is  used  for  REDOIT,  and  a  =  25, 

mt 

p  =  50  for  RAFSON  routine. 

If  SETUP  determines  that  an  iteration  is  necessary,  the  solution  of 
these  equations  is  obtained  through  an  iteration  method  using  an  alternating 
procedure  between  the  Newton-Rafson  method  and  binary  search  in  two 
dimensions. 

Twelve  iteration  steps  are  first  attempted  with  the  Newton-Rafson 
method;  the  initial  guess  on  a  and  8  is  either  specified  by  the  user  or  is 
determined  internally  if  the  user  so  desires.  If  this  method  fails  to  converge 
a  closer  approximation  is  obtained  for  the  initial  guess  using  one  step  of  the 
binary  search,  and  twelve  more  steps  are  performed  with  the  Rafson  method. 
If  this  still  fails,  another  closer  initial  guess  is  determined  with  two  steps  in 
the  binary  search,  and  twelve  more  Newton-Rafson  steps  are  taken.  Finally, 
^convergence  is  still  not  attained,  ten  steps  of  the  binary  search  are  taken. 

Subroutine  RAFSON 

This  subroutine  (see  Figure  33)  iterates  for  the  solution  to  the  equa¬ 
tions 

,  4(ct-a2  (q+8+3) _  . 

2  *  - 

(atl)  (8+1)  (a+g+4) 

3(a+B+3)  k  (a+B+2)2  +  fa+11  (B+lV(tf-3-4)l 
h  (a+l)(8t-l)(a+B+4)  (a+B+5) 

by  the  Newton-Rafson  method. 

An  initial  guess  is  made  by  the  SETUP  routine ,  with  this  (a,  8) » the 
first  partial  derivatives  of  these  equations  are  computed: 
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Figure  33.  Program  111:  Subroutine  RAFSON 


With  these  ,  the  increments  for  d.  and  g  are  determined  using  the 
relations  of  the  Newton-Rafson  method: 


'  ^(0,0-^ 

9^ 

W 

Sfu 

B2(a,0-32 

2 

33 

aa" 

93 

dP2 

aa 

*93 
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and 


P2  9  -  \  -  sP2/Sa  la 


Then  new  values  of  9  and  3  are  computed  by 
a=  a  +  Aa,  8  =  3  +  A3  » 

and  B  (a,  3) »  80  (a ,  8)  are  computed.  If  these  do  not  satisfy  the  criterion  of 
1  6 

0. 01  percent,  the  process  is  repeated  up  to  a  maximum  of  twelve  times  each 
time  it  is  called. 

Subroutine  REDOIT 

In  the  event  that  Newton-Rafson  iteration  fails  to  converge ,  this  sub¬ 
routine  (see  Figure  34)  is  used  to  compute  a  guess  for  RAFSON  that  will  be 
closer  than  the  last.  The  routine  is  also  used  to  attain  convergence  if  over¬ 
flow  or  divide  checks  makes  RAFSON  diverge. 

The  procedure  is  a  binary  search  in  two  dimensions.  An  initial  guess 
is  computed  by  SETUP.  Each  time  after  the  first,  the  last  value  is  used  as 
the  guess. 

The  initial  8  is  held  fixed,  and  a  is  incremented  by  35  percent,  until 
the  last  two  values  of  8  ^  (a ,  3)  fall  on  opposite  sides  of  3^  .  (K  after  20 
times  no  interval  containing  the  solution  exists ,  the  program  exits  and  skips 
this  case.) 

If  an  interval  is  found  successive,  linear  interpolations  are  performed 
until  the  computed  8^  is  within  0. 1  percent  of  the  desired  value  (maximum 
of  10  times).  This  completed,  the  same  thing  is  done  with  8  ,  holding  the 
newly  computed  a  fixed  until  82  *s  satisfied. 
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'  gr-'—" 


A  *  a 
B  >  0 
AX  ■  FIXED  o 
BX  •  FIXED  0 

AA  •  o  COMPUTED  FOR  STRADDLE 
BB  ■  0 

AAA  *  o  COMPUTED  FOR  INTERPOLATION 
BBB  *  0 

B1 A  .BlAA.Bi  AAA  >  0,  FOR  CORRESPONDING  e 
B2B,B2BB,62 BBB  -0{  *  0 


FBETAI  IS  A  FUNCTION 
|  THAT  COMPUTES  THE  EQUATION 
,  4  ( a-B  )8tafflS) 
{o+IH/J+lHa+/9+4r 


fit 


VINTRP  IS  THE 
LINEAR  INTER¬ 
POLATION  FUNCTION! 


6t20  1 

216 

INTERPOLATE 

COMPUTE  BIAAA 

FOR:  AAA 

FOR 

BETWEEN  AS  AA 

INTERPOLATION 

SET  A  OR  AA  EQUAL 
TO  AAA  SO  THAT  THE] 
TWO  BETAl'S 
STRADDLE  THE 
DESIRED/}, 


Figure  34.  Program  111:  Subroutine  REDOIT 
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With  the  new  (a ,  g) ,  g^  (a ,  0  is  computed.  If  the  results  are  within 
0. 1  percent  .of  ^  ,  set  the  indicator  to  and  return;  if  not,  the  process  is 
repeated.  The  number  of  times  the  process  is  repeated  is  a  function  of  the 
number  of  times  RAFSON  has  been  used; 


1  time 

2  times 

2  times 

4  times 

3  times 

8  times 

4  times 

12  times 

NOTE;  On  the  last  trial  for  REDOIT,  a  1  percent  criterion  instead  of  0.1 
percent  is  used.  If  convergence  is  not  attained  after  the  maximum  number  of 
steps,  set  the  indicator  to  1  and  return. 

The  routine  GASHUN  performs  the  computations  for  the  normalized 
beta  curve. 

Subroutine  GASHUN 

This  subroutine  (see  Figure  35)  computes  the  mean,  standard  deviation, 
and  tail  end  points  for  the  normalized  beta  curve.  The  mean 

q+l 

mN  oo-g+2  ’ 

and  standard  deviation 

a  _  VfcfrD  18+1) 

^  (a+fi+2)  Va+e+3 

are  computed. 

Then  the  following  integral 


is  to  be  evaluated; 
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A-13 


Figure  35.  Program  111:  Subroutine  GASHUN 


■100 


where 


bL-ma* 

bH-min  1^  +  40^ 


Simpson’s  rule  is  used  to  evaluate  the  area.  The  formula  for  the  area 
under  the  curve  over  a  particular  segment  is: 

A  =  |  h  [Y(l)  +  4Y(2)  +  2Y(3)  +  4Y(4)  +  Y(5)  ]  , 

where  the  Y's  are  five  Y-values  taken  at  five  points  within  this  segment  on  the 
X-axis. 

The  area  for  each  of  NINT  segments  (NINT  is  an  input)  is  computed  and 
added  to  the  total  of  all  previous  ones  to  get  the  total  area.  Then  the  seg¬ 
ments  are  accumulated  until  their  total  is  greater  than  A  times  the  total 
area  (A  is  an  input).  The  X  of  the  last  segment  is  the  lower  bound  (XNCL) 
for  the  tail  desired.  The  upper  bound  (XNCH)  is  determined  similarly. 
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